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ABSTRACT 

The method developed by RafFelt (1990a,b,c) to estimate a possible increase in the 
standard values of the helium-core mass at the tip of the red giant branch, M^, from 
properties of the color-magnitude diagrams of Galactic globular clusters is employed. 
In the present study, we revise and update Raffelt's database, including also constraints 
from RR Lyrae pulsation, and find that a small increase, by AM^ ~ 0.01 ± 0.015 Mq, 
cannot be ruled out with the present data and evolutionary models. Our new upper 
limits on AM^ are less restrictive than those previously obtained by Raffelt, as are 
the corresponding constraints on novel astroparticle phenomena which may affect the 
evolution of low-mass red giants. Within the estimated uncertainties, however, the 
standard values of may also be acceptable. Raffelt's method does not rule out a 
low envelope helium abundance in globular cluster giants, though again the standard 
values are compatible with the available constraints. The influence of a non- solar ratio 
for the a-capture elements upon these results is also investigated. In addition, we 
review several aspects of the input physics employed in red giant stellar evolutionary 
calculations, with the purpose of evaluating possible sources of uncertainty in the 
value of the helium-core mass at the helium flash that is obtained from evolutionary 
computations, such as: heat conduction by electrons in the degenerate core; Coulomb 
effects upon the Equation of State; triple-a reaction rates and screening factors; 
neutrino emission rates, both standard and enhanced by a possible non-zero magnetic 
moment; stellar rotation; microscopic element diffusion; and energy losses by axions 
and Weakly Interacting Massive Particles. 

Subject headings: Galaxy: Globular Clusters: General — Stars: Evolution — Stars: 
Giants — Stars: Horizontal-Branch — Stars: Interiors — Stars: Population II 
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1. Introduction 

The determination of the mass of the degenerate helium core of low-mass red giant stars at 
the tip of the red giant branch (RGB), M^, is a very important astrophysical problem for a variety 
of reasons. controls the luminosity of the star near the RGB tip, which is important for the 
determination of the rate of mass loss and thus the horizontal-branch (HB) morphology (e.g., Fusi 
Pecci & Renzini 1976). The luminosity of the RGB tip is important in distance-scale arguments 
as well (Lee, Freedman, & Madore 1993). Also, controls all the subsequent evolution on 
the HB and beyond, and is a crucial factor determining the sense of evolution of the RR Lyrae 
variables (e.g., Sweigart & Gross 1976): increased values would drive redward evolutionary 
paths for RR Lyrae variables in Oosterhoff type II globular clusters (GCs), and thus help account 
for the Sandage period-shift effect (Sandage 1982, 1990, 1993) and related problems in terms of 
the evolutionary scheme of Lee, Demarque, & Zinn (1990), as recently emphasized by Catelan 
(1992a,b, 1993a,b). Perhaps even more importantly, the luminosity of the HB depends markedly 
on Mc, and so does the distance (and hence age) scale that is most commonly employed for the 
GCs in the halo of the Galaxy and nearby systems. 

The evolutionary values obtained by Sweigart & Gross (1978, hereafter SG78) have been 
confirmed by several authors (e.g., Straniero & Chieffi 1991; VandenBerg 1992a,b). However, some 
computations have resulted in values which are not in clear agreement with the SG78 results (e.g., 
Harpaz & Kovetz 1988; Alongi et al. 1991; Raffelt & Dearborn 1987). In particular, Mazzitelli 
(1989, hereinafter M89) has claimed the SG78 values to be systematically too small by ~ 0.02 Mq, 
which he attributed to the inadequacy of the "shell-shifting" technique (Harm & Schwarzschild 
1966) for shifting the hydrogen-burning shell and advancing the chemical composition profile in 
an evolving RGB star and to an inappropriate choice of input physics by SG78. Very recently, 
Sweigart (1994) has analyzed in great detail the numerical aspects of the problem, and found that 
the Mc values do not hinge strongly on the shell- shifting algorithm employed. On the other hand, 
Sweigart suggests that the M89 prescriptions for implicitly advancing the chemical composition 
profile may lead to an overestimate of the helium-core mass, and thus should not be employed. 

Even though the numerical aspects of the problem of RGB computations have been put on 
a solid basis by Sweigart's (1994) analysis, not so much attention has been devoted to the input 
physics. In particular, Sweigart has not addressed the non-standard effects which may play a role 
in determining the core mass at the helium flash. 

It is the purpose of the present report to give an updated account of several aspects of the 
input physics that controls the core mass at the helium flash. We shall divide the present study 
into four main sections. In Sect. 2, Raffelt's (1990a,b,c) method to evaluate a possible increase 
in Mc is critically readdressed, on the basis of recent observational and theoretical results. The 
method suggests that a AM^ > cannot be ruled out with the present data, and may actually 
be favored by them. In Sect. 3, relevant standard ingredients of RGB evolutionary calculations 
are addressed: electron heat conductivities, the Equation of State (EOS), reaction rates, and 
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standard neutrino losses. In Sect. 4, the present status of important ingredients of non-standard 
evolutionary calculations, stellar rotation and microscopic diffusion, is critically discussed. In 
Sect. 5, exotic stellar energy loss mechanisms are considered, as caused by the presence of a 
non-zero neutrino magnetic moment, a "light" axion, and Weakly Interacting Massive Particles 
(WIMPs). Finally, our main conclusions are summarized in Sect. 6. 



2. AMc from the observations and stellar evolution: RafFelt's method revisited 

As already stated, the value of has a marked impact on the luminosities of the RGB tip 
stars and the RR Lyrae variables. Also, controls the lifetimes of the stars on the RGB and HB 
evolutionary phases and the RR Lyrae pulsation properties. 

Recently, Raffelt (1990a,b,c) has employed the available observational data and stellar 
evolution models to evaluate AMc = M^''^^''^ - M^^'''^'^'''"^. After an extensive analysis, Raffelt 
(1990a) arrived at two expressions for the increase in values [his eqs. (37) and (38)], and 
from these concluded that the increase in cannot be very large, with a likely 1 — a limit of 
AMc ~ 0.02 M0 [cf. his eq. (39)]. 

A revision of Raffelt's (1990a,b,c) analysis at the present time is appropriate for several 
reasons. First, the author has employed for the RR Lyrae variables a mean magnitude of 
(My) = 0.68 mag, which several new lines of evidence now suggest may be too faint (e.g., Saha et 
al. 1992; Catelan 1992a,b, 1993a,b; Simon & Clement 1993; Cacciari & Bruzzi 1993; Sandage 1993 
and references therein; Fernley 1994; see, however. Storm, Carney, & Jones 1994). In particular, 
the distance modulus of the Large Magellanic Cloud (LMC) determined from independent distance 
indicators, such as the analysis of the SN1987A circumstellar ring (e.g., Panagia et al. 1991; 
Crotts, Kunkel, & Heathcote 1995), gives strong support to this result (Walker 1992a). Some of 
these points have been recently addressed by Fernley (1994) and van den Bergh (1995). Second, 
RR Lyrae pulsation has not been addressed by Raffelt, even though it is known that the core mass 
influences the RR Lyrae properties in a very significant way (see, e.g., Catelan 1992a,b, 1993a,b 
and Castellani & De Santis 1994 for recent discussions). Finally, in Raffelt's original analysis 
several theoretical relations which are strictly applicable only to zero-age HB stars have been 
employed, while his sample of GCs includes several objects with extreme HB morphologies, where 
these variables may be well-evolved (Lee 1990; Catelan 1992a,b, 1993a,b; Caputo et al. 1993). 

We thus update Raffelt's (1990a,b,c) analysis in the present section, taking these problems 
into consideration. 
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2.1. Theoretical relations 

In RafFelt's (1990a,b,c) study, extensive use has been made of the relations derived by Buzzoni 
et al. (1983), which aimed at a semi-empirical determination of the helium abundance in GC stars 
through the so-called "iZ-method" (Iben 1968b). We shall similarly employ those relations in the 
present analysis. 

We assume that the helium flash does not induce mixing of metals and helium to the star's 
envelope (Iben & Renzini 1984; Renzini & Fusi Pecci 1988). This assumption has yet to be checked 
in detail by self-consistent evolutionary and hydrodynamic computations, which however are 
strongly dependent on initial and boundary conditions, as well as on the treatment of convective 
energy transport in the core (see, e.g., Deupree & Wallace 1987 for a discussion). 

2.1.1. The RR Lyrae magnitude level 

The dependence of the HB absolute visual magnitude at the RR Lyrae level on the envelope 
helium abundance Y, metallicity Z, and AM^, is given by the following relation: 

Y = 0.561 - 0.285 {M^^) + 0.046 log Z - 0.285 Arr - 2.083 AM^, (1) 

where Arr represents the deviation of the RR Lyrae luminosities from the zero-age predictions 
at logTeff = 3.85. A bolometric magnitude Mbd = 4.74 mag has been adopted for the sun. 
This follows from the analysis of recent satellite data, as summarized by Guenther et al. (1992), 
Sackman, Boothroyd, & Kraemer (1993), and Bahcall & Pinsonneault (1995), and the formula 
of Lang (1986) to convert luminosities into bolometric magnitudes. In Rafl"elt's (1990a) original 
expression [his eq. (9)], = 4.72 mag had been used instead. 

2.1.2. The brightness of RGB tip stars 

In order to estimate the absolute magnitudes of stars at the tip of the RGB, an assumption is 
required as to the distance moduli of the corresponding clusters (Frogel, Persson, & Cohen 1983). 

For our present purposes, we determine the magnitude difference between the HB at the RR 
Lyrae level and the RGB tip, AMy , and compare the observed values with the theoretical 
predictions (Rafl"elt 1990a). For the latter, following Buzzoni et al. (1983) and Rafl"elt (1990a), we 
flnd 

Y = 1.455 - 0.227 AM^'""^^ + 0.089 log Z - 0.227 Arr + 0.909 AM^. (2) 
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2.1.3. Star number counts and the R-method 



From the theoretical relations provided by Buzzoni et al. (1983) for the lifetimes of giants and 
HB stars, the envelope helium abundance may be estimated as follows: 



where R represents the ratio between the number of RGB stars brighter than the HB level and 
the total number of HB stars, including the RR Lyrae variables. 



The value of Mc controls the luminosity, mass, and temperature of a star inside the instability 
strip (e.g., Sweigart & Gross 1976), and thus, through the period-mean density relation, also its 
pulsation period, with an increase in generally leading to an increase in the period. 

Comparison between the standard models and the observations shows that the computed 
periods appear too small. Catelan (1992a,b, 1993a,b) has suggested that a systematic error in 
the evolutionary values may be responsible for this. However, the dependence of the periods 
on the effective temperature is large, which complicates the comparison between the models and 
the observations. On the observational front, small errors in the reddening can be of importance 
in "fixed-temperature" arguments (cf. Caputo & De Santis 1992 and references therein). On the 
theoretical side, the location of the blue and (especially) red edges of the instability strip is another 
delicate issue, and great theoretical effort is now being devoted toward their accurate evaluation 
(Bono & Stellingwerf 1994 and references therein). Since the RR Lyrae "mass-to-light ratio," 
A = logL — 0.81 logM, may be computed from the observational data in a reddening-independent 
way (Caputo & De Santis 1992), it may be more safely compared with the theoretical predictions 
at a fixed effective temperature for stars close to the zero-age HB. 

We have obtained the theoretical predictions for the zero-age HB mass-to-light ratio A at an 
effective temperature logTefr = 3.85 from the models of Sweigart, Renzini, & Tornambe (1987). 
The results displayed in their Table 4 may be fitted by the following equations: 



Y = 0.146 + 0.437 log R - 0.013 log Z - 0.144 Arr + 0.306 AM, 



(3) 



2.1.4. Pulsation properties of RR Lyrae stars 



logM 



1.122 - 0.082 log Z + 0.191 Yus + 1-294 M, 



(4) 



and 



logi 



0.451 - 0.026 log Z + 1.96 Yus + 3.197 M, 



(5) 



where M and L are the stellar mass and luminosity, respectively, both in solar units. Yms 
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represents the helium abundance prior to the main- sequence phase. With the aid of equations (1) 
and (6) from Raffelt 1990a, we finally obtain 

Y = -1.146 + 0.750 ^3.85 - 0.010 log Z - 1.548 AM^. (6) 

2.2. Observational relations 

2.2.1. The RR Lyrae magnitude level 

For the absolute magnitude of the HB at the (mean) RR Lyrae level, we employ Walker's 
(1992a) recent analysis of old GCs in the LMC with significant numbers of RR Lyrae variables. 
Assuming (m - M)o = 18.50 for the LMC (e.g., Panagia et al. 1991), Walker (1992a) finds 

(M^^) = 0.15 [Fe/H] + 0.73. (7) 

In the present work, we adopt the revised distance modulus for the LMC recently obtained by 
Crotts et al. (1995) through detailed analysis of the geometry of the SN1987A circumstellar 
envelope, (m — M)o = 18.58 ± 0.13, and convert the (My^) — [Fe/H] relation to the corresponding 
{M^^) - [Fe/H] one with the aid of Dor man's (1992) prescriptions (cf. Sandage 1993), so that 

{M^^) = 0.184 [Fe/H] + 0.68. (8) 

We adopt an uncertainty in values obtained from this equation of ±0.16 mag. 

It is important to note that, according to the data compiled in Walker's (1992b) Table 7, the 
Walker 1992a cluster sample is not biased in an evolutionary sense, since it does not appear to 
contain objects with extreme HB morphologies. Also, the adopted distance modulus for the LMC 
appears to be supported by several independent distance indicators, such as Cepheids and Miras, 
as discussed by Walker (1992a). 

2.2.2. The brightness of RGB tip stars 

When studying the brightness difference between the RGB tip and the HB at the RR Lyrae 
level, Raffelt (1990a), in his detailed analysis, has included some objects which have too few RR 
Lyrae variables and thus cannot be safely employed in a comparison with the zero-age HB relations 
within the instability strip. This is due to the possibility that, in clusters with too few RR Lyrae 
variables, these stars are evolved from the blue (or red) zero-age HB, and thus not representative 
of the zero-age HB at the RR Lyrae level (Lee et al. 1990; Lee 1990; Catelan 1992a,b, 1993a,b; 
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Caputo et al. 1993). Thus, we have repeated RafFelt's analysis excluding the clusters in his Table 2 
with HB morphologies outside the range -0.7 < {B - R)/{B + V + R) < +0.7, where B, R, and V 
are the number of stars to the blue, to the red, and within the instability strip, respectively (data 
taken from Lee, Demarque, & Zinn 1994). NGC 2808, which has the peculiarity of a dramatic 
internal second-parameter effect, has also been removed. We are left with the following objects: 
NGC 1261, NGC 1851, NGC 3201, NGC 5272, NGC 5904, NGC 6121, and NGC 7006. The 
metallicities have also been revised, following Lee et al. (1994). According to Raffelt's equations 
(22) to (25), and adopting a^hs = 0.1 mag, we then find 

(AM^^P"^^) = 4.09 ± 0.15, (9) 
for a mean metallicity ([Fe/H]) = -1.47 ± 0.15. 

2.2.3. Star number counts and the R-method 

Under the same restrictions as in the previous subsection, the data in Table 5 of Raffelt 
(1990a) associated to the number ratio counts have been reanalyzed. Clusters satisfying the HB 
morphology criterion were NGC 1851, NGC 3201, NGC 4147, NGC 5272, NGC 5904, and NGC 
6121. We find that, for such clusters, 

(log i?) = 0.123 ± 0.037, (10) 
for a mean metallicity ([Fe/H]) = -1.49 ± 0.17. 

2.2.4. Pulsation properties of RR Lyrae stars 

Employing a relation analogous to Caputo & De Santis's (1992) equation (12) (cf. Catelan 
1995) in connection with the observational data for M3 ([Fe/H] = —1.66) provided by Sandage 
(1990) in his Table 5, we find that 

A3.85 = 1.80 ±0.01. (11) 

In obtaining this value, we have employed Catelan's (1993a) procedure to search for the lower 
envelope of a distribution. In the present case, we assume that the zero-age HB corresponds to 
the lower envelope thus obtained. Sandage (1990) has argued that for this cluster this is probably 
a good approximation, as supported by synthetic HB models (cf. Catelan 1993a and references 
therein). Independent application of this method to the cluster NGC 3201 confirms this result. 
We caution, however, that the above error bar does not include systematic uncertainties, as might 
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be present if the temperature scale employed by Caputo & De Santis (1992) is not fully adequate. 
Further details will be given elsewhere (Catelan 1995). 



2.3. Theory and observations confronted 

To compare the theoretical and observational relations obtained above, it is necessary to 
evaluate the quantity Arr appearing in some of them. Raffelt (1990a,b,c) has neglected such 
terms in his analysis. In the present work, we follow the discussion by Catelan (1992a), where 
the variation of the HB magnitude dispersion at the RR Lyrae level as a function of [Fe/H] in 
Sandage's (1990) sample of GCs has been studied. 

Since 



AM^^ = 0.751 + 0.235 [Fe/H] (12) 
is Catelan's (1992a) prescription for the run of the HB thickness with cluster metallicity, and since 



(t?) = (41 ± 15)% (13) 

is the mean "concentration parameter" (measuring the mean deviation of the variables from 
the lower envelope of the distribution) for those clusters studied by him with intermediate HB 
morphology, it follows that 



Arr ~ 0.31 + 0.10 [Fe/H]. (14) 

We now insert these results into the theoretical relations given in Sect. 2.1. For the RR Lyrae 
magnitudes, for consistency with the observational data for R and AMy , we adopt the value 
given in equation (8) for a metallicity [Fe/H] = —1.48 ± 0.16, and find 

Y = (0.250 ± 0.049) - 2.083 AM^, (for M^^) (15) 

where (t{M^^) = 0.16 mag and (7(Arr) = 0.06 mag have been assumed. In a similar vein, from 
the number ratio counts (iZ-method), we find 



F = (0.218 ± 0.018) + 0.306 AMc. (fori?) (16) 
From the magnitude difference between the HB and the RGB tip, one finds 
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Y = (0.207 ± 0.039) + 0.909 AM^. 



(for AM\ 



rtip-RR 



) 



(17) 



Finally, the mass-to-light parameter A yields 



Y = (0.238 ± 0.008) - 1.548 AMc 



(for A3.85) 



(18) 



Each of these four equations defines a band in the Y — AM^ plane, as plotted in Fig. 1. As 
can be seen, all four bands intersect at a relatively narrow "compatibility region" in the middle of 
the diagram, around (F, AM^/Mq) ~ (0.221 ± 0.017, +0.011 ± 0.014), which defines the preferred 
solution for the envelope helium abundance and the increase in M^. Regions depicted in lighter 
gray tones in this diagram define progressively less likely solutions. 

Analysis of this diagram reveals that, even though the data are generally compatible with 
the standard values, an increase by ~ 0.010 — 0.015 Mq (or possibly larger) cannot be ruled 
out, and in fact may appear to be favored. It is interesting to note that the narrow band running 
from upper left to lower right in this diagram, coming from equation (18), is of great importance 
in constraining the allowed solutions. Such a constraint, obtained from the analysis of RR Lyrae 
pulsation properties, had not been previously taken into consideration in the ignition argument 
(Raffelt 1990a,b,c). Systematic uncertainties were not included in the analysis of A3. §5, though, 
so that the actual "compatibility region" may be slightly wider. In any event, our 1 — a limit 
for AMc is appreciably larger than Raffelt's, posing less stringent constraints on non-standard 
physical phenomena that may be affecting the evolution of giant stars (cf. Sect. 5). 

Lower envelope helium abundances are also consistent with Fig. 1. As already mentioned, 
a low Y is favored by the evolutionary interpretation of the Sandage period-shift effect. A more 
detailed analysis of the helium abundance problem can be found in Catelan 1992b, where the 
difficulties in accounting for an abnormally low Y in GCs were discussed. In particular, a low 
helium abundance may not appear compatible with primordial "^He nucleosynthesis arguments, 
both theoretical and observational (e.g., Pagel et al. 1992; Pagel & Kazlauskas 1992; Olive, 
Steigman, & Walker 1991; Olive & Steigman 1994). Fig. 1 may also suggest that a slightly larger 
Y would imply a AMc closer to zero. Furthermore, relaxation of the pulsation bounds would allow 
considerably higher Y values. At any rate, it must be emphasized that a low envelope helium 
abundance in an evolved star does not necessarily imply a low initial abundance for this star, 
since diffusive mechanisms may decrease the envelope abundance along a star's lifetime. 



2.4. 



Influence of an overabundance of the a-elements 



The above estimates have been made assuming solar-scaled abundances for the GC stars, 
i.e., [a/Fe] = 0. However, it is well known (see Wheeler, Sneden, & Truran 1989 for a review) 



- 10 - 



that halo stars present an overabundance of the a-elements, typically by a factor ~ 3 (i.e., [a/Fe] 
~ 0.48). Therefore, the above analysis may require modifications due to this specific effect. 

Until now, the only extensive study dealing with the evolutionary implications of an 
overabundance of the a-elements in low-mass stars is the one by Chieffi, Straniero, & Safaris 
(1991) and Salaris, Chieffi, & Straniero (1993). According to Chieffi et al., the a-enhanced models 
may be obtained from the solar-scaled models on the basis of a straightforward rescaling of the 
latter in terms of the overall metal abundance Z, following a relation of the type 

log Z = [Fe/H] - 1.7 + log(0.579 / + 0.421), (19) 

where / = lO^"/^'']. For the evolutionary stages that precede the HB phase, these results have 
been confirmed by Chaboyer, Sarajedini, & Demarque (1992). 

Following the Chieffi et al. (1991) prescriptions, we have repeated the above analysis 
assuming [a/Fe] = 0.48. Our results are summarized in Fig. 2. Even though the impact of the 
a-enhancement upon the four individual relations between Y and AM^ is quite significant — 
larger, in fact, than suggested by the Raffelt & Weiss (1992) analysis — the overall picture does 
not appear to change appreciably. In particular, the "compatibility region" now appears slightly 
narrower, being found around {Y, AM^/Mq) ~ (0.223 ± 0.012, +0.008 ± 0.011). 

Although the agreement between the model predictions and all the available observational 
constraints does seem remarkably satisfactory, we emphasize that an increase in reaching 
perhaps as much as 0.015 Mq cannot be ruled out. Since Sweigart (1994) has very convincingly 
argued against an increase in due to numerical problems in state-of-the-art evolutionary 
computations, it should be worthwhile to check the input physics in more detail. 

3. Standard evolutionary effects, standard physics 

In the present section, we shall analyze important ingredients of standard evolutionary 
calculations, aiming at an evaluation of their present uncertainties and the corresponding impact 
upon AMc- 

3.1. Electron thermal conductivities in the interior of RGB stars 

As one advances inward in the interior of an RGB star, electron conduction is progressively 
more efficient as an energy transport mechanism in comparison with radiative transfer, and 
becomes the dominant mechanism in the electron-degenerate helium core (see, e.g., the discussion 
accompanying Fig. 1 in Iben 1968a). 

It is common practice, in stellar evolution calculations, to employ either the conductivities 
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that have been presented by Hubbard & Lampe (1969, hereinafter HL69) — as in SG78 — or 
the more recent results by Itoh et al. (1983, hereinafter I'^MSS) — as in M89. To our knowledge, 
however, a discussion of the hypotheses and uncertainties involved in the adoption of either 
formalism has not been presented in the literature. We suspect that this is due to the fact that, 
since the work by Iben (1968a), an understanding has grown that "electron conduction opacities 
in degenerate cores are so small anyway that the core is isothermal in any case" (quotation from 
Raffelt & Weiss 1992). However, Sweigart et al. (1987) have shown that, by multiplying conductive 
opacities by a factor = 0.5, may increase by 0.0087 — 0.01 Mq (cf. their Table 6). Since 
this is in the ballpark of the "compatibility regions" in Figs. 1 and 2, some exploration of this 
topic should be of interest. 

"Conductive opacity," k^, is usually defined, in terms of the more familiar conductivity Ac, as 
follows: 



4acT3 , ^ 

= W 

where a is the radiation density constant, c is the speed of light, p is the density, and T is the 
temperature. The usefulness of such a definition is that it places the heat conduction and radiative 
transfer equations in identical forms (e.g., Clayton 1983). 

In Fig. 3, a comparison is made of the conductive opacities available in the literature for the 
physical conditions prevailing in the center of RGB stars. For illustrative purposes, we plot the 
conductive opacities along the RGB evolution of the M = 0.9 Mq, Yms = 0.20, Z = 10""* model 
of SG78, based upon the log and log values displayed in their Table 2. For convenience, we 
have assumed that = 1. 

In Fig. 3a, the conductive opacities are directly displayed. The classical analysis of HL69 
is represented both by interpolative equation (A4) of Iben 1975 (hereafter Ib75) and by direct 
linear interpolation on their published tables. Equation (A28) of Ib75 represents Canuto's (1970) 
relativistic results. An interpolative scheme has been suggested by Ib75 [his eqs. (A31) and (A32)], 
whereby the HL69 results are valid for log p < 6, and Canuto's for log p > 6.3. For the Yakovlev 
& Urpin (1980) case, we have adopted equation (19) of Blinnikov & Dunina-Barkovskaya (1992). 
Equation (7) of I3M83 has been employed as representative of their analysis, and the quantum 
corrections suggested by Mitake, Ichimaru, & Itoh (1984) were evaluated on the basis of their 
equation (8). Finally, the contribution to conductive opacity due to electron-electron interactions, 
as calculated by Urpin & Yakovlev (1980), was obtained from equation (27) of Blinnikov & 
Dunina-Barkovskaya (1992), with the aid of equation (10) by Timmes (1992). 

In Fig. 3b, the ratio of conductive opacities from the different authors to the Ib75 
approximation to the HL69 results is shown. Several interesting features are apparent. 



First, the Canuto (1970) formulation coincides closely with the Ib75 approximation to the 
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HL69 results for p ~ 10^ g cm~'^, so that the relativistic corrections proposed by this author 
are of no practical importance for the (essentially classical) core. As a consequence, the curve 
labeled "Hubbard & Lampe" matches almost perfectly the one obtained from the Ib75 proposed 
interpolation between Canuto's (1970) and HL69's results all across the diagram of Fig. 3a, 
which is the reason why the latter is not displayed. Similarly, the quantum corrections that have 
been obtained by Mitake et al. (1984) are also very small. Third, the Yakovlev & Urpin (1980) 
opacities are somewhat larger than those due to HL69. Last but not least, it is also noted that 
the electron-electron contribution to the conductive opacity, as prescribed by Urpin & Yakovlev 
(1980), may be quite significant near the RGB tip. 

Conductive opacities, as displayed in Fig. 3, never disagree by more than ~ 45%. Caution is 
needed, however, before concluding that present uncertainties in of order 0.01 Mq from this 
source are an upper limit. In particular, it should be noted that: 



The formulation of I^M83 and Mitake et al. (1984) has been developed aiming at applications 
to the more extreme conditions characterizing white dwarfs and neutron stars. It is not 
expected to he especially accurate for the interior of an RGB star. In Fig. 3b, a curve has 
been drawn to show the values achieved by the parameter 

Z^e^ 

" (21) 



r-ksT 



where e is the electron charge, ks is the Boltzmann constant, Z is the atomic number, and 
r is the ion-sphere radius, r = (3/47m;)^/'^, with n; the ion number density. T measures the 
strength of the electrostatic interaction among the plasma ions. [A related quantity. A, is 
also displayed, and will be discussed in the next section.] For T ^ 178, the ions form a period 
structure (i.e., a crystal), whereas for F ^ 2 a Boltzmann gas is still present. We see that F 
never exceeds 0.81 in the RGB core. The calculations of PM83 and Mitake et al. (1984) '^'''^ 
valid only for 2 < F < 160, i.e., for the liquid phase, as explicitly stated by these authors. 
Thus, one should keep in mind that employing their results for low-mass RGB stars actually 
involves an extrapolation to the Boltzmann gas conditions that better characterize the ions 
in such stars; 

To shed light on the uncertainties involved in such extrapolations, Canuto's (1970) relativistic 
results were extrapolated towards log p < 6, as shown in Fig. 3. The agreement with the 
IIL69 opacities clearly deteriorates for the lower densities; 

The Itoh et al. (1984) results should not be employed in RGB models, since they refer to 
even more extreme conditions, 171 < F < 5000 (there is also a recent erratum to this paper 
which should be referred to). The recent analyses by Itoh & Kohyama (1993) and Itoh, 
Hayashi, & Kohyama (1993) again refer exclusively to the crystalline lattice phase, and the 
present authors are unaware of an extension of this work to the regions of interest for RGB 
calculations (see also Itoh 1992); 
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• In their compilation of input physics for white dwarf models, Blinnikov & Dunina- 
Barkovskaya (1992) have pointed out that the calculations of I'^M83 and Yakovlev & Urpin 
(1980) do not include the effects of electron-electron interactions on the conductivities. As 
already noted, the e — e contribution to the opacity for the relevant conditions may be quite 
significant (cf. Fig. 3). This point may have been neglected in some recent evolutionary 
calculations. The e — e and e — i contributions to the opacity add linearly, in a first 
approximation; 

• Nandkumar & Pethick (1984) have also carried out an analysis similar to the one by I'^M83, 
giving results for T > 4. They find that inclusion of physical ingredients neglected by the 
latter authors, related to the influence of the electrons upon the effective interactions among 
the ions and between ions and electrons, can lead to differences of a few percent in the 
conductivities at low densities. Inspection of their Fig. 1, however, discloses that the effect 
referred to is negligible at ~ 10^ g cm~'^; 

• The Ib75 approximation to the IIL69 conductive opacities appears to underestimate the 
values obtained by direct interpolation on IIL69's Table 2 in the regions of interest for RGB 
interior computations. Approximation formulas of this kind are sometimes employed in stellar 
evolution codes. Our comparison shows that, whenever possible, accurate interpolation 
routines should be employed to obtain the opacities from the original references. 

From these considerations, a few additional remarks are needed. One of the major differences 
between the M89 and SG78 analyses is that in the former the more recent (I'^M83) conductive 
opacities were adopted. If the extrapolation carried out is appropriate, then Fig. 3 shows the 
conductive opacities of I'^M83 to be somewhat smaller than those obtained by direct interpolation 
on IIL69's tables (but larger than the values obtained from Ib75's proposed approximation!). 
This being the case, the core mass values attained with the I'^M83 opacities should be slightly 
larger (i.e., by a few times 0.001 Mq), since it would be easier for the helium core to cool and 
the helium flash would take place later. Part of this difference might be removed if the e — e 
interaction were properly taken into account. Note, however, that the Urpin & Yakovlev (1980) 
e — e curves in Fig. 3 are themselves extrapolations, since these authors have provided results 
only for p > 10^ g cm~'^. One run of Mazzitelli's models with the IIL69 conductivities has indeed 
resulted in a smaller M^, although it is difficult to isolate the impact of upon this result since 
several other ingredients were also changed in the same run. Thus, the disagreement between M89 
and SG78 may be partially ascribed to the particular formulations for conductivities employed by 
these authors. 

It appears clear to us that useful work could still be developed in the field of conductive 
opacities with the purpose of estimating more precisely the value of M^. In particular, evaluation 
of the inaccuracies involved in extrapolating the I'^M83 calculations toward P < 1 is needed. While 
this is not done, we are inclined to suggest that preference should be given to the IIL69 results in 
population II evolutionary calculations. 
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To close this section, we note, after the discussion by Iben (1968a), that inaccuracies in 
radiative opacities near the core boundary may also influence values in an important way. 
This should also be kept in mind when updating the input physics in an evolution code. 



3.2. Coulomb effects upon the EOS 

Coulomb efl"ects upon the EOS were not included in the SG78 computations. According to 
Harpaz & Kovetz (1988), Coulomb interactions lead to a reduction in by ~ 0.005 Mq, as a 
consequence of a reduction in the pressure leading to higher core temperatures, and thus to a flash 
that takes place earlier. The efl"ect is not very large because at the conditions prevailing in RGB 
cores the EOS is dominated by electron degeneracy. 

Coulomb corrections, when implemented in stellar evolution codes, are generally based on 
the Debye & Hiickel (1923) formalism or some modiflcation thereof (see, for instance, Straniero 
1988 for a discussion and references). However, according to Rogers (1994), the Debye-Hiickel 
correction to the EOS is only valid for 



" -(wF^ « 

Outside this range, more detailed treatments are required. As shown in Fig. 3b, A > 1.1 in an 
RGB core, so that the Debye-Hiickel formalism breaks down and more sophisticated analyses are 
warranted to accurately model the evolution of RGB properties. 

Only very recently (Chaboyer & Kim 1995) has the detailed OPAL EOS (cf. Rogers 1994 
and references therein) been incorporated into evolutionary computations. However, these authors 
could not evolve their models up the RGB due to numerical difficulties. Extending their analysis 
to the RGB phase should be of great importance for the determination of the actual influence of 
Coulomb efl"ects upon the values attained at the flash. ^ 



^In this regard, it is important to note that the determination of fully self-consistent AT^ ages with the OPAL 
EOS must await the computation of detailed RGB sequences, since HB evolutionary characteristics — luminosities 
in particular — depend on Afc, as discussed in Sect. 2. As an example, a decrease in Afc by Si 0.005 Mq due to EOS 
effects (Harpaz &; Kovetz 1988) would increase h.V ages by ~ 3% — 4%, since it would make the HB appear fainter. 
The age reduction in Chaboyer &; Kim 1995 (by ~ 6% — 7%) is a semi- empirical effect, i.e., it is strictly valid only if a 
stellar evolution-independent calibration of the HB magnitudes is employed, or if the turnoff magnitude is estimated 
from some other means. 
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3.3. Reaction rates 

The helium flash in the core of an RGB star occurs when the triple-a process, in which three 
"^He nuclei are transformed into a ^^C nucleus, is ignited, with the building of a trace equilibrium 
abundance of ®Be nuclei (e.g., Clayton 1983; Castellani 1985; SG78). 

Near T ~ 10® K, the rate of the triple-a reaction may be given as 

eaa oc T^o (23) 

[see, e.g., eq. (5-105) in Clayton 1983]. On the other hand, Sweigart (1994) has recently estimated 
that the maximum temperature within the degenerate helium core of an RGB star near the RGB 
tip increases at a rate dlog Tmax/dMgh ~ 1.4 Mq^, due to the increasing mass of hydrogen that 
is converted to helium at the hydrogen-burning shell. Thus, in order to delay the helium flash 
in such a way as to enable an increase in the core mass by, say, 0.015 Mq, an increase in the 
maximum core temperature by AlogT ~ 0.021 should occur prior to helium ignition. This would 
imply that the present reaction rates are too high by A log £30, ~ 40 X 0.021 = 0.84, i.e., a factor of 
IqO.sa ^ 7_ [Sweigart's (1994) estimate of an error by a factor of 16 corresponds to an increase in 
the core mass by ~ 0.02 Mq.] Thus, as pointed out by Sweigart (1994), quite a large error in the 
reaction rates would have to be present to cause an important increase in M^- 

The reaction rates most commonly employed for the two-step triple-a process are those given 
by Fowler, Caughlan, & Zimmerman (1975), Harris et al. (1983), Caughlan et al. (1985), and 
Caughlan & Fowler (1988). The latter two references give identical rates for this reaction, while 
some difl"erences exist with respect to the others. It may be of some interest to reinvestigate the 
origin of such difl"erences, since in the classical study by SG78 the earlier rates provided by Fowler 
et al. (1975) were used for the triple-a reactions. 

We have computed the ratio between the more recent rates for the triple-a process and those 
provided by Fowler et al. (1975). For the conditions near the helium flash (T ~ 10® K), we flnd 
that 

e3«(1988) ~ 1.22 x e3„(1975), (24) 

which implies a decrease in the core mass by only ~ 0.002 Mq. The stated uncertainty in either 
reaction rate is somewhat larger than the difl"erence between the two, albeit much smaller than, 
say, a factor of two. [To be sure, there is a second branch in the tabulated reaction rates which 
is uncertain by a large factor, but this becomes more important than the main branch only for 
T ^ 51 X 10® K, being thus irrelevant for the flash conditions.] 

Oberhummer et al. (1994) have recently analyzed the triple-a process from a difl"erent 
perspective. These authors have shown that even extremely small variations in the strength of 
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the effective nucleon-nucleon interaction (i.e., by ~ 0.1% — 0.2%) may have a large impact upon 
the reaction rates. The authors interpret their results in terms of the change that is introduced in 
the Boltzmann factors exp(— £^R/kBT) due to the shifts in the resonance energies of the ®Be and 
^^C nuclei in the two steps of the reaction. For instance, a change in the resonance energy of the 
dominating O2 state for the ®Be(a,7)^^C reaction from the present value of E^i = 0.2875 MeV 
to = 0.2517 MeV would, by itself, change the reaction rate by a factor of 60 (cf. their Table 
1). We estimate, from equation (4-195) of Clayton (1983), that a change by a factor of 7 in the 
reaction rate would require a smaller resonance energy for this state by ~ 6%. Though this may 
seem small, we believe it is unlikely that the experimental determinations could be wrong by such 
an amount. 

Finally, it is important to recall, after Tarbell & Rood (1975), that the so-called screening 
factors of the nuclear reactions are another relevant ingredient that influences the value of the core 
mass at the RGB tip, since they can considerably modify the unscreened nuclear reaction rates. 

In modern stellar evolution computations, the screening factors usually employed are those 
by DeWitt, Graboske, & Cooper (1973) and Graboske et al. (1973), who have presented a general 
expression that is claimed to be valid in the so-called weak, intermediate, and strong screening 
regimes. More recently, Yakovlev & Shalybkov (1989) have readdressed in detail the problem of 
ion and electron screening in dense stellar interiors. These authors claim that their analysis is 
more accurate than the one by DeWitt et al. (1973) and Graboske et al. (1973), due to their use of 
more reliable calculations of the relevant Coulomb thermodynamic quantities. To our knowledge, 
the Yakovlev- Shalybkov prescriptions have not been incorporated into modern evolutionary 
computations, so that it would be worthwhile to compare these more recent results with those 
from Graboske et al. 1973. 

For the present calculations, we assume, for simplicity, p = 10^ g cm"^, T = 10® K, and 
Yc = 1. For the gaseous, strongly electron-degenerate case, ion and electron screening for the 
triple-a process may be obtained from equations (184) and (186) in Yakovlev & Shalybkov 1989, 
respectively. For an enhancement factor given by 

/ = exp(^), (25) 
we flnd, assuming intermediate screening to be operating. 

Hi = 1.396, /i = 4.039, (26) 
where the subscript i denotes ion screening. As to electron screening, we flnd 



= 0.042, /e = 1.043. 



(27) 
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Given that /e is very close to unity, we conclude that electron screening is weak. 

Ion screening may also be obtained from equation (19) of Graboske et al. 1973, with the aid 
of their Table 4. Assuming = — which amounts to considering extreme degeneracy and no 
electron screening — the usual assumption in RGB interior calculations (cf. SG78), we find 

H; = 1.326, /i = 3.766. (28) 

[Note that the enhancement factor obtained from eq. (A6) in Tarbell & Rood 1975 for intermediate 
screening is ~ 2 times larger than this result, even though their expression for weak screening, eq. 
(A5), agrees perfectly well with the value obtained from the Graboske et al. 1973 expression.] 

This result is very similar to the one obtained from the Yakovlev & Shalybkov (1989) 
formalism, and we conclude that the differences between the Graboske et al. (1973) and Yakovlev 
& Shalybkov analyses are of minor relevance for the quantitative evaluation of M^. 

Unfortunately, the treatment of screening in the non- degenerate conditions that characterize 
the H-burning shell employed by Sweigart & Gross (SG78 and references therein) was not described 
in detail in their series of papers. Checking the influence of this ingredient in future calculations 
should also be of some interest. 

3.4. Standard neutrino losses 

Uncertainties in the standard fitting formulas for the neutrino emission rates are another 
potential source of error in the determination of Mc, as already emphasized by Chieffi & Straniero 
(1989). These authors note, for instance, that the standard energy loss rates due to photo, pair, 
and plasma processes, as given in the fitting formulas by Munakata, Kohyama, & Itoh (1985), 
are valid for the range < log p < 14, 8 < log T < 10 whereas in RGB interiors, where typically 
log ~ 6, energy losses by neutrinos begin to be a relevant cooling mechanism well below 
log T = 8. Note that evolutionary computations are commonly based upon the rates periodically 
improved by the Japanese group (cf. Itoh et al. 1992 and references therein). 

Only very recently, though, a fitting formula has appeared (Haft, Raffelt, & Weiss 1994) 
that is sufficiently reliable for application in stellar evolution calculations aiming at a precise 
evaluation of (cf. Altherr & Salati 1994). Comparison of the new plasma emission rates 
with those previously obtained by several different groups reveals surprisingly large differences 
(cf. Figure 8 in Haft et al. 1994). Haft et al. have concluded that the correct emission rates in 
the parameter regime relevant for the RGB in GCs are larger by about 10% — 20% than those 
of (former) stellar evolution calculations, implying an increase in the core mass by ~ 0.005 Mq. 
The SG78 calculations, in particular, neglect neutral- current effects (Weinberg 1967; Salam 1968), 
which may have led to an underestimate of their values (Ramadurai 1976; SG78). 
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4. Non-standard evolutionary effects, standard physics: stellar rotation and 

microscopic diffusion 

On the basis of the Mengel & Gross (1976) rotational models, Renzini (1977) has obtained the 
following relation between the core mass increase AM^ and the initial rate of stellar rotation lij: 

AMc = 1.25 X 10^a;2-i6^ (29) 

where is given in rads~^. Thus, a stellar rotation rate as small as ~ 2.2 X 10~'*rads~^ might be 
sufficient to produce an increase in as large as 0.015 Mq. 

The physical argument behind the Mengel & Gross (1976) results is quite simple: in the 
rotating stellar interior of an RGB star, besides degeneracy pressure working against further 
compression of the stellar material, there is also a centrifugal component working against the 
compression of the core and the eventual ignition of the triple-a reactions, thus retarding the 
helium flash and leading to an increase in over the non-rotating model predictions. 

Following these theoretical arguments and based upon equation (29), much interpretative 
work on the properties of rotating RGB and (ensuing) HB structures has been developed (e.g., 
Renzini 1977; Fusi Pecci & Renzini 1976, 1978; Castellani, Ponte, & Tornambe 1980; Castellani 
& Tornambe 1981; Sweigart et al. 1987; Renzini & Fusi Pecci 1988; Shi 1995). Peterson (1983, 
1985a,b) provided observational evidence that evolved GC stars do rotate while on the HB, and 
also suggested a very interesting trend of increasing rotation rate for clusters with bluer HB 
morphologies, which appeared to go in favor of a signiflcant influence of rotation upon stellar 
evolution on the RGB and beyond. 

Recently, more sophisticated stellar evolution codes that take into account the efl"ects of 
rotation upon the evolving structure were developed. In these models, the influence of angular 
momentum loss from the stellar surface and of angular momentum redistribution within the 
stellar interior were consistently taken into account (Deliyannis, Demarque, & Pinsonneault 1989; 
Sofla, Pinsonneault, & Deliyannis 1990; Pinsonneault, Kawaler, & Demarque 1990; Pinsonneault, 
Deliyannis, & Demarque 1991a). This represented a clear improvement over the earlier models by 
Mengel & Gross (1976), where a simple one- dimensional approximation with angular momentum 
conservation in individual mass shells was adopted. 

The predictions of these recent stellar evolution models were very difl"erent from those in the 
simple approximation by Mengel & Gross (1976). Deliyannis et al. (1989), in particular, in their 
Figs. 3 and 4, compare the amount of angular momentum remaining in their models near the 
zero-age main sequence and turnofl" regions with the corresponding predictions in the Mengel & 
Gross case, concluding that the latter (but not the former) are inconsistent with the observed 
surface rotation velocities of HB stars and the observational limits for rotation near the turnofl". 
Most importantly for our present purposes, they point out that their models with difl"erential 
rotation are unable to preserve sufficient angular momentum in the core as to cause a signiflcant 
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increase in over the values expected for non-rotating models, even though the core retains 
sufficient angular momentum that the observed HB velocities can be accounted for. They claim, as 
a consequence, that rotation cannot affect the core mass and the HB morphology in a significant 
way, and even speculate that the correlation obtained by Peterson (1985b) between stellar rotation 
of HB stars and HB morphology could be driven by metallicity and age effects. 

It appears likely, however, that the issue has not been definitely resolved yet. The amount 
of initial angular momentum in population II stars, the angular momentum loss law and several 
physical instabilities that may be responsible for angular momentum redistribution in the stellar 
interior, are some of the aspects of the problem that may be improved upon. Besides, there are 
many free parameters in these models that must be adjusted in a semi-empirical way. Some of the 
most salient uncertainties have been pointed out by Michaud & Charbonneau (1991, especially 
their Sects. 6.6 and 7.5). We note as well that these models have not been evolved beyond the 
main-sequence turnoff region, so that the conclusions about more advanced phases have yet to be 
established in a more accurate way. For instance, the CNO abundances in RGB stars of GCs do 
not appear to be accounted for by these models without treating mixing processes that should 
take place on the RGB phase itself (Pinsonneault et al. 1991b). On the other hand, Sofia et al. 
(1990) have argued that since the rotational models do not appear to be in conffict with any 
particular observational data (including the light element abundances), any additional physical 
ingredient to be added should be of sufficiently small relevance. They also point out that many of 
the model results are not particularly sensitive to their choice of free parameters, this appearing 
to be the case also for the core mass at the flash (cf. Deliyannis et al. 1989). 

With respect to the suggestion that age is a key variable responsible for the differences in 
rotational velocities among HB stars in GCs, there now appears to be at least one example that 
goes against this idea. While it is clear that M3 and M13 have very different stellar rotational 
velocities (Peterson 1985b) and HB morphologies for nearly identical metallicities, it has recently 
been argued that these two clusters do not have appreciably different ages (Catelan & de Freitas 
Pacheco 1995).^ 

It is well known that element diffusion may be responsible for generating sufficiently 
high molecular weight gradients as to inhibit rotational instabilities that provoke mixing (e.g., 
meridional circulation). In this case, the good agreement between rotational model predictions and 
the observations, in what concerns the light element abundances, might be rather fortuitous. In 
the case of low-mass stars, diffusion is mainly caused by gravitational settling, given that thermal 
diffusion is small for the conditions prevailing in such stars (cf. Deliyannis, Demarque, & Kawaler 
1990), as is radiative levitation (Michaud, Fontaine, & Beaudet 1984). 

Much work has recently been developed to investigate the impact of microscopic diffusion on 



^ After this paper had been completed, the authors became acquainted with the latest investigation by Peterson, 
Rood, &; Crocker (1995) on the correlation among stellar rotation, HB morphology, and abundances on the RGB. 
Their analysis of the clusters M13 and M3 appear to confirm our arguments. 
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stellar evolution predictions for population II stars. This has been motivated by the possibility of 
obtaining significant reductions in estimated GC turnoff ages (Stringfellow et al. 1983) and (in 
diffusion and rotation studies alike) by the need to estimate the cosmological lithium abundance 
from the observed Spite & Spite (1982) Li plateau for metal-poor dwarfs (Deliyannis et al. 1989, 
1991). 

However, apart from the discussion in Proffitt & VandenBerg (1991) regarding the impact 
of He diffusion upon the luminosities of HB stars and hence the ages that are estimated on 
the basis of the AV method, the influence of diffusion upon has not been studied in detail. 
Chaboyer et al. (1992), in their recent analysis of GC ages under helium diffusion, assume that the 
standard results remain valid in the diffusive approximation, provided a lower initial Y is adopted. 
This assumption appears to be valid (see below) only insofar as the envelope helium abundance 
decreases uniformly as a function of [Fe/H] due to the influence of helium diffusion. Preliminary 
evidence in favor of this assumption has been presented by Proffitt & VandenBerg (1991), but 
more calculations for more metal-poor compositions are needed to corroborate their results. 

Proffitt & VandenBerg (1991) argue that "the core mass that is ultimately attained at the 
helium flash is mainly a function of the envelope helium content," and also that "the evolution of a 
giant which has a decreased envelope helium abundance due to diffusion should be essentially the 
same as that which did not undergo gravitational settling but which simply had a lower Y when 
it was formed." [Recall that the timescale for diffusion is such that the amount of He diffused 
is appreciable in relatively long-lived evolutionary phases only, but not on the RGB itself (cf. 
Deliyannis et al. 1990 for a detailed discussion).] 

Assuming this argument to be correct, it follows that the change in due to microscopic 
diffusion can be obtained from the relation 



^ = -O.24M0, (30) 

as prescribed by the SG78 and M89 models. From Fig. 3 in Proffitt & VandenBerg 1991, 
l^standard " l^diffusion G [0.006, 0.024], implying a AMJMq G [+0.00015, +0.006] for [Fe/H] = 
— 1.26 and [0/Fe] = 0.55. Shi (1995) has recently claimed that "Mc will be slightly increased by 
0.005 Mq with respect to that of standard models due to the helium settlement into the core." 
This is likely to be an upper limit, particularly in higher metallicity stars with thicker convection 
zones, since diffusion in convection zones is strongly inhibited. 

Note that the Proffitt & VandenBerg (1991) argument is quite general, and may be extended 
to whatever processes that may lead to a change in the envelope helium abundance. For instance, 
meridional circulation, driven by rotational instabilities, may increase the envelope helium 
abundance of an RGB star and lead to a smaller M^. Deliyannis et al. (1989) dismiss this (subtle) 
effect. Deep mixing, often invoked to explain abundance anomalies in RGB stars (see Freeman & 
Norris 1981 for a review), may play a similar role if started early enough on the RGB. 
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The interplay among these various mixing mechanisms is far from being well established, and 
it is possible indeed that mutual inhibitions take place. The effects of helium diffusion can inhibit 
rotational instabilities, which can inhibit diffusion, which in turn may be much less effective 
if overshooting or other type of convective instability is taking place (but see also Proffitt & 
Michaud 1991). Deliyannis & Demarque (1991) have concluded, on the basis of a comparison of 
their Li isochrones with the observed high-temperature end of the Li plateau, that their models 
overestimate diffusion in population II stars (but see also Proffitt & Michaud 1991 and Chaboyer 
et al. 1991), so that inhibiting mechanisms must be effective; Swenson (1995) has recently pointed 
out that mass loss on the main-sequence phase may improve the agreement, even though it is 
quite unclear whether the required rates, ~ lO~^^M0yr~^, are attained in such stars. Much more 
theoretical and observational work is clearly necessary. Very recently, Chaboyer et al. (1995) have 
attempted to include rotation, diffusion, and overshooting simultaneously in their solar models, 
and found the final results to be significantly different from those achieved when these different 
processes are individually considered. When comparing their predictions of solar oscillation modes 
with the observations, they find indications of important missing ingredients in their treatment of 
rotation. Similar calculations should be carried out for population II stars. 

One may find interesting that Figs. 1 and 2 suggest a simultaneous increase in and 
decrease mY — in the same sense as expected on the basis of equation (30) above. However, to 
account for an increase in by 0.015 Mq, equation (30) would suggest a decrease in Y by 0.06 
— which is probably unrealistic. 

To close the present section on the role of standard physical ingredients in non-standard 
evolutionary computations, we mention in passing that the recent results by Castellani & 
Castellani (1993) represent strong evidence against a significant influence of mass loss on the RGB 
upon Mc- 

5. Non-standard evolutionary effects, non-standard physics 

In this section, some representative novel physical ingredients are briefly discussed in the 
framework of RGB stellar evolution: a non-zero magnetic moment for the neutrino, the axion, and 
Weakly Interacting Massive Particles. Even though there are important reasons to regard these 
three exotic mechanisms as the main candidates for playing actual roles in the evolution of stars, 
several other novel phenomena have been proposed whose impact on stellar evolution may be of 
relevance, such as millicharged particles (Davidson & Peskin 1994 and references therein), a "fifth 
force" (e.g., Grifols, Masso, & Peris 1989), etc., but which shall not be addressed here. The reader 
is referred to Raffelt 1990b for a critical discussion of some of these possibilities. 
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5.1. A non-zero neutrino magnetic moment 

One proposed explanation for the "solar neutrino problem," motivated by a possible 
anticorrelation between solar magnetic activity and the measured solar neutrino flux on Earth (cf. 
Oakley et al. 1994 and references therein), is that electron neutrinos originating from the solar 
interior change their nature when traversing the solar convective zone. More speciflcally, it has 
been suggested that, due to an anomalously large magnetic moment, left-handed neutrinos would 
turn into right-handed neutrinos when reaching the magnetic fleld in the solar convective zone, 
and thus escape detection when reaching the Earth. Estimates by Voloshin, Vysotskii, & Okun 
(1986) suggest that a magnetic moment as high as /x^.^ ~ 10~^°/xb would be needed (where /xb is 
Bohr's magneton), although more detailed considerations about the magnetic flelds in the solar 
interior may allow smaller values for /x^.^ (Oakley et al. 1994). 

Presently, laboratory experiments (cf. Vogel & Engel 1989 and references therein) appear 
unable to constrain the size of /x^.^ to better than an order of magnitude above this value. On 
the other hand, astrophysical arguments may be able to improve this bound. For instance, it 
has been claimed (e.g., Fukugita & Yazaki 1987) that the cosmological "^He abundance constrains 
the magnetic moment to /x^.^ < 5 X 10~^^/xb. Similarly, it has been argued that the right-handed 
neutrinos escaping SN1987A could be reflipped again to left-handed neutrinos, and then be 
detected on Earth: analysis of the neutrino signal from SN1987A, as monitored by the Kamiokande 
detectors, puts stringent constraints on the size of the neutrino magnetic moment (Notzold 1988 
and references therein), which however are mo del- dependent (see, for instance, Blinnikov & Okun 
1988 and Vogel & Engel 1989 for discussions and additional references). Rafl"elt (1990b) provides 
an excellent review of the subject. 

Another important astrophysical implication of a non-zero neutrino magnetic moment would 
be an increased rate of cooling of stellar degenerate cores. This is due to the increased (over 
the standard rates discussed in Sect. 3.4) energy losses from the dominant decay of plasmons 
Tpl VeVe (Bernstein, Ruderman, & Feinberg 1963; Sutherland et al. 1976). 

This aspect of the problem has been most recently addressed by Rafl"elt (1990a), Rafl"elt 
& Weiss (1992), and Castellani & Degl'Innocenti (1993). In particular, in the last two papers 
the efl"ect of increased cooling due to the increased plasmon decay rate has been incorporated 
into detailed stellar evolution codes. Rafl"elt & Weiss (1992) compare their results with those by 
Castellani & Degl'Innocenti (1993), flnding good agreement to within the approximations adopted 
in the latter study. 

Both Rafl"elt & Weiss (1992) and Castellani & Degl'Innocenti (1993) point out that relations 
of the form depicted in Sect. 2 are in fact only flrst-order approximations to the real case. In fact, 
the relation between core mass and RGB luminosity actually depends sensitively on the type of 
cooling involved: the standard relation between these two quantities is not really universal. This 
implies that observational bounds on a possible core-mass excess will translate into slightly weaker 
bounds on /x^.^. Hence, we assume that the non-standard increase in from Sect. 2 translates 
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into bounds on /x^.^ which are less restrictive by ~ 10"^Vb (RafFelt & Weiss 1992). 

Employing RafFelt's (1990a) equation (45), corrected by RafFelt & Weiss (1992) to read [cf. 
their eq. (18)] 

AMc = 0.013 ^^/XB, (31) 

we find that a non-standard increase in by 0.015 Mq could be accounted for by a neutrino 
magnetic moment of 



/X,, ^ 1.15 X IO-^Vb- (32) 

This value applies directly to the diagonal magnetic or electric dipole moments of Dirac neutrinos, 
or to the transition moments of Majorana neutrinos. For the transition moments of Dirac 
neutrinos, such a value should be multiplied by a factor 2~^/^ (RafFelt 1990b). 

It is important to note that there appears to be no stronger semi-empirical constraint 
on the magnetic moment of the neutrino. In particular, Wang (1992a) and Blinnikov & 
Dunina-Barkovskaya (1992, 1993) have recently argued that the cooling of hot white dwarfs places 
the most reliable constraints on the neutrino magnetic moment, /x^.^ < 4 X 10~^^ /xb (Wang 1992a), 
which would enable an increase in by as much as 0.05 Mq. This is ruled out by Figs. 1 and 2 
above. Notzold's (1988) limit on fj,^^ from SN1987A considerations, fj,^^ < (1.5 ± 0.5) X 10"^^ /xb, 
imply AMc < (0.020 ± 0.007) Mq. 

Barbieri & Fiorentini (1987a,b) have recently argued, from simple theoretical considerations, 
that the neutrino magnetic moment should be smaller than fj,^^ ~ 10"^'^ /xb. If this is true, neutrino 
magnetic moments produce negligible effects upon the He-core mass (i.e., AMc ^ 10"'^ Mq). 



5.2. Axions 

The axion is a hypothetical, (almost) massless particle, which appears as a pseudo-scalar 
Goldstone boson in a class of theoretical particle models. The axion gives an elegant solution to 
the problem of CP conservation in strong interactions, and is a cold dark matter candidate (cf. 
RafFelt 1990b for a comprehensive review). 

As discussed by RafFelt & Dearborn (1987), astrophysical considerations stand out as a 
unique opportunity of constraining the properties of the physical axion, since it may have a strong 
influence upon important evolutionary stages of stars, depending however on the details of the 
specific axion model. 

Hadronic axions couple to ordinary nucleons but not to electrons, and arguments related 
to the solar age and solar helium abundance (Dearborn, Schramm, & Steigman 1986; RafFelt & 
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Dearborn 1987) constrain its mass to be 

m^^lleVR. (33) 

In a large class of grand unified theories, R = 1. 

As in Sect. 5.1, we must check whether the standard core mass-luminosity relation is valid in 
the case of cooling of RGB interiors by axions. It is interesting to note that, contrary to the case 
in which only cooling by neutrinos is considered, cooling by axions produces a central helium flash 
(cf. Fig. 7e in Rafl"elt & Dearborn 1987). 

From analysis of Table III in Rafl"elt & Dearborn 1987, we flnd that, for a 1.3 Mq star, the 
axion mass which leads to an increase in by 0.069 Mq (ma — 2.1 eV) leads to an increase in 
the RGB tip luminosity by Alogi ~ 0.34. A similar increase in M^, according to equation (10) in 
Rafl"elt 1990a, would lead to a brighter RGB tip by Alogi ~ 0.324. Since we are concerned with 
smaller core mass increments than the one considered in this speciflc example, we may neglect 
entirely this small difl"erence: the direct efl"ect of a possible axion cooling upon the RGB tip 
luminosity is therefore not relevant. 

A very difl"erent efl"ect, however, concerns the duration of the helium-burning phase. When 
we discussed neutrino cooling, we assumed that its sole implication upon helium-burning stars 
would be through its influence upon M^. This is a sufficiently accurate approximation, since 
neutrino cooling is very inefficient for the physical conditions prevailing in the interior of an HB 
star. Energy losses by hadronic axions, on the other hand, dramatically increase while the star is 
on the HB (e.g., Raffelt 1990b), so that there is a further, direct influence upon the duration of 
the HB phase. 

The following comparison illustrates the importance of this effect. The axionic evolutionary 
models by Raffelt & Dearborn (1987) show that an axion mass of — 2.1 eV leads to an 
increase in the core mass by 0.069 Mq, and to a reduction in the duration of the HB phase by 
AlogfHB — 0.88. Equation (16) in Raffelt 1990a, on the other hand, indicates that the reduction 
due to the quoted increase in alone should be by only AlogfnB — 0.20. The difference is due 
to the fact that, in the latter relation, the direct impact of axions upon the HB phase has been 
neglected. This point has also been addressed by Haxton & Lee (1991). 

We thus conclude that the relations obtained by Raffelt (1990a,b,c) involving the iZ-method 
arguments cannot be plainly employed to constrain a non-standard increase in (see also the 
discussion in Raffelt 1990b). However, the remaining relations employed in Sect. 2.1. can still 
be regarded as reasonable approximations to the core mass excess that might arise from axion 
emission on the RGB phase, which however should be checked against the HB lifetime constraint 
(see also Raffelt & Dearborn 1988). 

Interpolating on Table III of Raffelt & Dearborn 1987, one may roughly estimate the required 
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axion mass to produce a AM^ = 0.015 Mq to be 



TTia ~ 0.6 eV, (hadronic axions) 



(34) 



which is more than one order of magnitude smaller than the upper bound inferred from the solar 
age argument (see above). From Table IV in Raffelt & Dearborn (1987), however, this is found to 
imply a reduction in the duration of the HB phase by ~ 63%, which is rejected on the basis of star 
count arguments (e.g., Raffelt & Dearborn 1988). Therefore, hadronic axions cannot be entirely 



In addition to the hadronic axions, DFSZ axions (Dine, Fischler, & Srednicki 1981; Zhinitsky 
1980) are believed to play an important role in particle physics. In the case of the DFSZ axions, 
which do couple to electrons, the emission occurs mainly through Compton and bremsstrahlung 
processes, while the Primakoff effect is the dominant emission process in the case of hadronic 
axions (Raffelt 1990b). The latest results by Raffelt & Weiss (1995) on the evolution of a low-mass, 
low-metallicity RGB star suggest that a core mass increase by 0.015 Mq could be accounted for 
by DFSZ axions of mass 



Interestingly, the rate of increase in the pulsation period of the ZZ Ceti variable G117-B15A 
(Kepler et al. 1991), which appears too large in comparison with the expectations from CO white 
dwarf models, has motivated an accelerated cooling scenario (Isern, Hernanz, & Garcia-Berro 
1992) in which an excess cooling by axions of an estimated mass ma — (8 ± 3)/ cos^ /3 meV might 
be taking place. However, recent results by Vauclair (1994) suggest that the variability is not 
as large as previously reported, and also that no obvious periodicity is present. Therefore, we 
believe that this argument offers no compelling evidence for axion emission, although this case is 
illustrative of the kind of observations that could lead to the positive discovery of non-standard 
cooling. It should be noted, in addition, that the axion emission rates employed by Isern et al. 
may have been underestimated (Altherr, Petitgirard, & del Rio Gaztelurrutia 1994), so that their 
TTia value may have been overestimated accordingly. 

The mass value for the axion which would be needed to produce a AMc ~ 0.015 Mq is not 
excluded by analysis of the white dwarf luminosity function (Raffelt 1986; Raffelt & Weiss 1995; 
Altherr et al. 1994) nor by the limits imposed by type I supernova explosions (Wang 1992b). 
SN1987A considerations may not rule out ma values lying close to that found in equation (35) 
above either (Altherr et al. 1994). 



responsible for an increase in by as much as ~ 0.015 Mq. 



ma K, 7.4 X 10"^/ cos^ /3 eV. (DFSZ axions) 



(35) 
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5.3. Weakly Interacting Massive Particles 

Another proposed candidate for the dark matter in the Universe are the so-called WIMPs. A 
wide variety of such non-baryonic, hypothetical particles has been suggested in the literature. A 
special class of WIMPs, termed "cosmions," has been proposed with properties which would help 
account for several astrophysical problems simultaneously: the solar neutrino problem; the age 
discrepancy between the oldest GCs and the measured expansion rate of our Universe; the dark 
matter problem; and, more recently, the "hump" in the GC stellar luminosity functions in the 
subgiant region (see, e.g., Faulkner 1991 for a review). 

Since the existence of cosmions is presently only speculative, wide debate exists in the 
literature concerning their impact on stellar evolution. Different views on the problem have been 
published by Renzini (1987), Spergel & Faulkner (1988), Faulkner & Swenson (1988, 1993), Rood 
& Renzini (1990), Dearborn et al. (1990), Finzi (1991), Riazi (1991), Bouquet & Salati (1989), 
Bouquet, Kaplan, & Martin de Volnay (1989), Salati & Bouquet (1990), Martin de Volnay, 
Bouquet, & Kaplan (1990), Raffelt (1990d), etc. In general, these authors have computed the 
abundance of WIMPs in the stellar interior for different stellar evolutionary phases, assuming 
those particles to have been accreted from the dark Galactic component along the stellar lifetime. 

We are mainly interested in the impact of efficient energy transport by WIMPs upon the 
low-mass red giants. Bouquet & Salati (1989) and Bouquet et al. (1989) have tackled in detail 
the problem of cosmion accretion and evaporation from stellar structures, and predicted that 
in the equilibrium condition cosmions will be so extremely concentrated in the RGB core, as a 
consequence of its much higher density in comparison with the more external layers, that no 
impact is expected upon RGB properties. Rood & Renzini (1990), on the other hand, have pointed 
out that a non-negligible influence upon the core evolution may be expected, given the non-local 
nature of WIMP energy transport. Since WIMPs are efficient energy carriers, a delay in the 
helium ignition may indeed be expected — one should recall that effects of order 1% in are of 
relevance for our present purposes. Actual RGB evolutionary computations might be of interest 
in this regard. [Interestingly, Finzi's (1991) models suggest that a small core of hydrogen should 
develop at the very center of a WIMP-accreting RGB star, so that explosive pycnonuclear burning 
of hydrogen (and not thermonuclear burning of helium) would trigger the core flash, which would 
thus occur earlier and at a smaller value than in standard calculations.] 



* Actually, the problem of cosmion evaporation from the stellar interior is complicated by the fact that the relevant 
WIMP cross sections are either purely axial or purely scalar (or vector) (cf. Faulkner &; Swenson 1993). In the case 
that the cross-sections are purely axial (or spin- dependent), interaction with the stellar material would depend on 
the presence of protons, given that a-particles are spinless and other nuclei with spin are much less abundant. In 
this case, transport effects by cosmions would decrease as hydrogen were converted to helium in the stellar interior, 
and thus little impact would be expected upon Mc values. On the other hand, for scalar or vector interactions, 
transport effects would increase as the core grew richer in helium, so that a more significant impact on Mc would 
appear reasonable. 
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The impact of WIMPs on the evolution of HB stars has been critically addressed by Renzini 
(1987), Spergel & Faulkner (1988), Dearborn et al. (1990), and Raffelt (1990d). In the latter two 
papers, evidence is found that the HB configurations become unstable in the presence of WIMPs, 
and that thermal pulses take place on short timescales until near core helium exhaustion. But 
little influence upon the HB lifetime has been found, in contrast with the simple arguments by 
Renzini (1987). Implications are also expected for RR Lyrae pulsation and the HB luminosity, 
as discussed by Dearborn et al. (1990) and Rafl"elt (1990d). These authors have found that the 
HB becomes brighter by ~ 0.3 mag with respect to non-WIMPy models; however, the HB mass 
and metallicity explored by them are not very representative of GC stars, and evidence exists 
that the situation would be much more complex for stars that actually fall within the RR Lyrae 
instability strip. Note that the increase in HB brightness found by the above authors ensues 
irrespective of an accompanying increase in M^, so that the relations obtained in Sect. 2 again 
cannot be straightforwardly employed when constraining the WIMP parameter space. At any rate, 
Spergel & Faulkner (1988) have noted that these calculations depend on assumptions about the 
evaporation of WIMPs during the helium flash, when the core expands in a very short timescale, 
and the central density drops from ~ lO^gcm"'^ (on the RGB tip) down to ~ lO'^gcm"'^ (on 
the HB). Most of the authors who have considered the efl"ect of WIMPs on the HB evolution have 
assumed that WIMPs do not evaporate in the flash. 

Faulkner & Swenson (1993) have recently reviewed the evidence that the luminosity functions 
of well-observed GCs at the subgiant region cannot be reconciled with the predictions from 
standard stellar evolution. This topic has also been addressed quite recently by Stetson (1991), 
VandenBerg (1993), and VandenBerg & Stetson (1991). Interestingly, the theoretical luminosity 
functions of VandenBerg appear to actually favor the development of isothermal cores in GC stars 
as an explanation for the observed luminosity functions; Fig. 8 in Stetson 1991 is particularly 
impressive in this regard, showing that the "WIMPy" models lead to clearly improved matches 
with the observations in comparison with standard (as well as diflusive) evolutionary models. On 
the other hand, VandenBerg & Stetson (1991) have found that the "hook" feature in the CMD of 
the old open cluster M67 may not be accounted for by their WIMPy models, a discovery that has 
been termed by Faulkner (1991) as possibly "the most stringent stellar evolutionary test for the 
WIMP hypothesis." The hook feature is well accounted for by standard models. 

In spite of the remarkable theoretical efl"orts that have been made by the above workers, it 
should be pointed out that the likelihood of cosmions as a solution to the solar neutrino problem 
(and thus to other problems in astrophysics) has decreased considerably in the recent years, as 
more data appeared on the solar neutrinos from the several independent groups now carrying out 
experiments in this fleld. Bahcall (1994) has reviewed the data for the ^Be and ®B neutrinos and 
argued very strongly against the "astrophysical solution" to the solar neutrino problem (i.e., a 
decreased central temperature leading to a smaller neutrino flux from the sun, as expected in the 
WIMP hypothesis). Schramm (1994) has also studied the problem in great detail, and even though 
his statistical analysis does not deflnitely rule out the astrophysical solution to the problem, it 
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clearly favours explanations involving new neutrino physics — e.g., neutrino oscillations. Hata, 
Bludman, & Langacker (1994 and references therein) have strongly argued that the astrophysical 
solution is ruled out by comparison of the Kamiokande and Homestake solar neutrino data. 
The review by Berezinsky (1994) provides a useful entry to the recent literature, and in his 
Sect. 5 shows the trends expected in future, much more sensitive, solar neutrino experiments 
(SNO, Superkamiokande, Borexino) if the solar neutrino problem can be accounted for by the 
astrophysical solution — an idea against which Berezinsky as well strongly argues. 

Several authors have also noted (e.g., Krauss 1990; Turner 1993) that the expected annihilation 
cross sections lead to a WIMP abundance that is simply too low (both within the sun and in 
the Universe as a whole) for the most "interesting" cosmion candidates, representing another 
extremely serious challenge to the cosmion hypothesis. In fact, recent developments do not even 
offer compelling evidence for the presence of significant amounts of non-baryonic dark matter in 
the scale of galactic halos, as reviewed by Turner (1993) and Schramm (1994, particularly his 
Fig. 5). In reality, this may be viewed as an argument against a large impact of axions upon 
stellar evolution in the Galaxy as well. 

Irrespective of these difficulties to the WIMP hypothesis, one must recall that an independent 
satisfactory explanation to the anomalies in the luminosity functions of GCs has not been proposed 
as yet. Even if it turns out that WIMPs cannot be responsible for this anomaly, some other form 
of efficient central energy transport may have to be identified, if Faulkner & Swenson's (1993) 
arguments turn out to be correct. Implications for the core evolution of giants should then be 
studied in detail, checking whether an increase in over standard values is among the implied 
effects. 

6. Conclusions 

Revision and extension of Raffelt's (1990a,b,c) method and database to estimate 
a possible increase in helium-core mass values at the tip of the RGB over standard 
prescriptions from observational data and evolutionary models for RGB and HB stars 
suggest a AMc ~ 0.01 ± 0.015 Mq. Our 1 — a limits on AM^ are higher than those previously 
obtained, and the ensuing constraints on novel physical mechanisms correspondingly weaker. 
Perhaps somewhat surprisingly, the method suggests that such an increase should be accompanied 
by a slight decrease in the envelope helium abundance. Similar results obtain for a-enhanced 
compositions, assuming the Chieffi et al. (1991) scheme for treating the problem to be equally valid 
for RGB and HB stars. A high helium abundance for GCs, as recently advocated by Shi (1995), 
is ruled out by this method. A high Y would also be in conflict with the explanation that Lee 
et al. (1990) have advanced for the Sandage period-shift effect (Sandage 1982, 1990, 1993), since 
RR Lyrae variables in Oosterhoff type II clusters would not be well-evolved objects, as required in 
this scenario. It would also produce pronounced blue loops in the evolutionary tracks (Sweigart 
& Gross 1976), which lead to "bifurcations" in synthetic HB models (e.g., Catelan 1992b, 1993a) 
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that are not observed in real clusters. These important points have been missed in Shi's (1995) 
paper. The combination of a decreased Y (Lee et al. 1990, Lee 1990) and an increased (Catelan 
1992a,b, 1993a,b), on the other hand, would render the Yale evolutionary interpretation of the 
Sandage effect more attractive, although an explanation should then be found for this apparent 
anomaly in the helium abundance. A low primordial Y seems to be ruled out, both on theoretical 
and observational grounds (e.g., Pagel et al. 1992; Pagel & Kazlauskas 1992; Olive et al. 1991; 
Olive & Steigman 1994). 

Several aspects of the input physics involved in evolutionary calculations which may bear a 
direct influence upon values have been reviewed in the present article, with the purpose of 
identifying likely causes for a possible core mass increase. In particular, the following aspects have 
been discussed: conduction properties of dense stellar matter; Coulomb efl"ects upon the EOS; 
reaction rates; cooling by neutrinos and the influence of its (possible) non-zero magnetic moment; 
stellar rotation and diflusion; energy losses by hadronic and DFSZ axions and by WIMPs. 

Our analysis reveals that some of these aspects have already been sufficiently elaborated 
in recent work that an increase in by ~ 0.005 Mq is already indicated, as in the case of 
neutrino cooling (Haft et al. 1994). Efficient cooling by exotic particles (or some other mechanism) 
should also be considered seriously, after the challenging results of Faulkner & Swenson (1993) 
on the luminosity function of well-studied GCs: in particular, standard models appear incapable 
of accounting for the observed number of stars in the sub-giant region of the color-magnitude 
diagram. Uncertainties in the reaction rates appear small enough that changes in due to 
them may safely be regarded small. ^ Screening in the interior of an RGB star is also a secondary 
source of uncertainty; however, we were unable to ascertain the influence of screening upon the 
hydrogen-burning reactions taking place in the vicinity of the helium core. Although recent results 
tend to place rotation as an irrelevant ingredient insofar as values are concerned (Deliyannis et 
al. 1989), several sources of uncertainty still exist which would recommend caution before regarding 
this the flnal word on the subject. In particular, it is to be noted that the "age-driven" rotation 
effect on HB stars that has been speculated by these authors as an explanation to Peterson's 
(1983, 1985a,b) observed correlation between stellar rotation rate and HB morphology must be 
discarded if M13 (blue HB morphology and fast HB rotators) and M3 (even HB morphology and 
slower HB rotators) are really essentially coeval, as recently claimed by Catelan & de Freitas 
Pacheco (1995). Any effect connected to rotation or diffusion and leading to a change in surface 
helium abundances is also expected to lead to a change in the core mass value at the RGB tip; 
there is still much uncertainty in this fleld, where useful work could still be carried out. We 
point out, in addition, that one of the most commonly employed sources of conductive opacities 
for stellar interior calculations, I'^M83, actually refers to the liquid metal phase (P > 2) only, 
whereas in the RGB core one hardly achieves P > 0.8. Thus, an uncertain extrapolation from a 



^The surprisingly different reaction rates recently reported by Fowler (1993) for the triple-a process are probably 
due to his consideration of a very different physical environment. 
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particular physical phase (liquid) to a different one (gas) is usually carried out. It is also worth 
noting that I'^M83 have neglected the e — e contribution to the conductive opacities (Blinnikov & 
Dunina-Barkovskaya 1992), an inadequate approximation for the RGB interior. To the present 
authors, the results of HL69 would appear more reliable for the calculation of population II RGB 
sequences, but care should be taken to employ sufficiently reliable interpolation schemes: we find 
the Ib75 approximation formulas to underestimate the HL69 opacities in the regions of interest 
for RGB computations. It would be of great interest to extend the results by the Japanese group 
toward smaller V. Employing the I'^M83 expressions in their present form to RGB interiors, a 
larger core mass is expected in comparison with those based upon IIL69, which may help explain 
the disagreement between the SG78 and M89 results. Finally, we call attention to the fact that 
the Debye-Hiickel correction to the EOS is not valid for the dense conditions characterizing RGB 
cores, so that a more refined treatment is required for the computation of accurate values. 

The authors are indebted to Drs. Georg Raffelt, Allen Sweigart, Ben Dorman, Don 
VandenBerg, and an anonymous referee for carefully reading drafts of this paper and the 
manuscript and for providing many important comments and/or suggestions. M. C. acknowledges 
financial support by FAPESP (grant 92/2747-8). J. E. H. has been partially supported by CNPq. 
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Fig. 1. — Graphical method to estimate Y and AM^ from the models and the observations. Four 
bands are displayed, obtained from analysis of the quantities My^, AMy^~^^, R, and A3. §5 as 
described in the text. There is a relatively narrow "compatibility region" in the middle of the 
diagram, around (Y, AMc/Mq) ~ (0.221, +0.011), where all of the four bands intersect, which 
thus defines the preferred solutions for Y and AM^. 



Fig. 2. — As in Fig. 1, except that an overall enhancement of the a-elements, [a/Fe] = 0.48, has 
been assumed. Note that the "compatibility region" is around (F, AMc/Mq) ^ (0.223, +0.008). 



Fig. 3. — Comparison among several different prescriptions for the conductive opacities along the 
evolution of the core properties of an RGB star. In panel (a), the conductive opacities are given, 
having been obtained from the indicated studies. In panel (b), the opacities are compared with 
the reference values from the Ib75 approximation to the IIL69 results. The curves are plotted for a 
star with evolutionary parameters {M/Mq,Yms, = (0.9, 0.20, 10""*). The core densities (given 
in cgs units) and temperatures were obtained from SG78. We have assumed a pure-helium plasma. 
The thick solid lines display the evolution of T, a non-dimensional parameter that measures the 
strength of electrostatic interactions in the medium (F ^ 2 for a gas, and F ^ 178 for a crystal); 
and A, which measures the range of validity (A ^ 0.2) of the Debye-Hiickel approximation to the 
Coulomb effects upon the EOS. 
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